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Prnp -1 ' mice lack the prion protein PrP c and are resistant to prion infections, but variable 
phenotypes have been reported in Prnp~l~ mice and the physiological function of PrP c remains 
poorly understood. Here we examined a cell-autonomous phenotype, inhibition of macrophage 
phagocytosis of apoptotic cells, previously reported in Prnp~ ! ~ mice. Using formal genetic, 
genomic, and immunological analyses, we found that the regulation of phagocytosis previously 
ascribed to PrP c is instead controlled by a linked locus encoding the signal regulatory protein a 
[Sirpa). These findings indicate that control of phagocytosis was previously misattributed to the 
prion protein and illustrate the requirement for stringent approaches to eliminate confounding 
effects of flanking genes in studies modeling human disease in gene-targeted mice. The pleth- 
ora of seemingly unrelated functions attributed to PrP c suggests that additional phenotypes 
reported in Prnp~l~ mice may actually relate to Sirpa or other genetic confounders. 



The cellular prion protein PrP c , encoded by 
the Prnp gene, is tethered to the membrane of 
most mammalian cells by a glycosylphosphati- 
dylinositol anchor. Conversion and aggregation 
of PrP c into a misfolded conformer (termed 
PrP Sc ) triggers transmissible spongiform encepha- 
lopathies, also termed prion diseases (Aguzzi 
and Calella, 2009). Disparate functions have 
been ascribed to PrP c on the basis of pheno- 
types described in Pmp~'~ mice (Steele et al., 
2007; Linden et al., 2008), yet none of these 
functions has been clarified mechanistically, and 
their validity was frequently challenged. 

All currendy available Pmp~'~ lines were gen- 
erated using embryonic stem (ES) cells derived 
from the 129 strain of Mus mnsculus. Typically, 
chimeric founder mice were crossed withWT 
(Prnp 1 """") mice of the C57BL/6 strain (B6;Table 1). 
The resulting mice carrying the targeted allele 
were maintained as a mixed B6/129 back- 
ground (B6129) or backcrossed to inbred strains 
(mainly B6), often for >10 generations, to gen- 
erate congenicB6.129 mice (the period in B6. 129 
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denotes congenesis, as opposed to mixed B6129 
mice). Although the latter approach progres- 
sively reduces the contribution of the 129- 
derived genome, Pn;p-linked genomic intervals 
may co-segregate with Prnp (Sparkes et al., 1986). 
Consequently, congenic Prnp'"""" and Prnp~'~ 
mice may differ at additional polymorphic loci 
(Smithies and Maeda, 1995; Gerlai, 1996). We 
hypothesized that co-segregation of linked 
genes may have confounded the attribution 
of functions to PrP c based on phenotypes ob- 
served in Pmp~'~ mice (Collinge et al., 1994; 
Lledo et al, 1996;Walz et al., 1999; Rangel et al, 
2007; Lauren et al., 2009; Calella et al., 2010; 
Gimbel et al, 2010; Ratte et al., 2011; Striebel 
et al.,2013). 

Here we selected a cell-autonomous pheno- 
type previously reported in congenic B6.129- 
Pmp ZM,ZlchI mice (Biieler et al., 1992) inhibition 
of phagocytosis of apoptotic cells (de Almeida 
et al., 2005) .We used RNA sequencing to identify 
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Table 1 . Prnp KO mouse lines analyzed in this study 


Prnp KO mouse line 


ES cells 


Origin of ES cells 


Strain of partner of 


Location of colony 


Reference 








chimeric mouse 






PfflpZrchl/Zrchl 


AB1 


129S7/SvEvBrd 


B6 


Zurich, Switzerland 


Bueleretal. (1992) 


PmpNgsk/Ngsk 


J1 


129S4/SvJae 


B6 


Nagasaki, Japan 


Sakaguchi et al. (1995) 


p rn pEdbg/Edbg 


E14 


129/Ola 


129/Ola 


Edinburgh, Scotland, UK 


Manson et al. (1994) 


p rnp GFP/GFP 


HM-1 


129/Ola 


B6 


Cambridge, MA 


Heikenwalder et al. (2008) 


p rn pRkn/Rkn 


E14 


129P2/0laHsd 


B6 


Wako-shi, Japan 


Yokoyama et al. (2001) 


p rn pZrchll/Zrchll 


E14.1 


129/OlaHsd 


B6 


Zurich, Switzerland 


Rossi etal. (2001) 


PfflpRcmO/RcmO 


HM-1 


129/Ola 


129/Ola 


Edinburgh, Scotland, UK 


Moore et al. (1995) 



The present study is based on the analysis of mice carrying seven independently generated Prnp-null alleles. Prnff db9/ab9 and Pmp"™ 0 /"™ 0 were always crossed to isogenic 
129/Ola mice, whereas all other Prnp~t~ mice were crossed to B6 mice and then kept on a mixed B6 and 129 background or further backcrossed to B6 or other strains. 



genes linked to Prnp and expressed in macrophages that may 
influence this phenotype.We report genetic and functional evi- 
dence that the regulation of phagocytosis previously ascribed 
to Pmp~ l ~ is instead controlled by the closely linked gene signal 
regulatory protein a (Sirpa; Matozaki et al., 2009). 

RESULTS 

RNA sequencing indicates Sirpa as a candidate 
modulator of phagocytosis 

Because hyperphagocytosis was reported in p rn p z 'M/zrM 
mice after extensive backcrossing to B6, we reasoned that any 
genetic confounders leading to this phenotype in B6.129- 
p rl] pZrdii/zniii 1TUCe WO uld need to be physically linked to Prnp, 
polymorphic between 129 and B6 strains, and expressed in 
macrophages. To identify such genes, we performed RNA se- 
quencing on B6A29-Prnp ZrM/Znh! and B6.129-Pr/7p"" / "" peri- 
toneal macrophages (pM$s). Of 11,586 genes, 305 were 
differentially expressed between the two groups (Fig. 1 A and 
Table SI). We then analyzed the reads for the presence of 
deviations from the B6 reference transcriptome (Fig. 1 B and 
Table 2) and applied a step-wise filtering strategy to narrow 
down potential confounding candidate genes (Fig. 1 C). 
B6.129-Prnp z ' chI/ZlTl ' 1 pM<J>s showed significantly more devia- 
tions from the B6 reference than did B6.129-Prnp"" /u " (Fig. 1 B 
and Table 2). By plotting the density of sequence variations 
in Bb.\29 -Prnp z " hllz,M pM$ mRNA against their physical 
location, we observed clustering of variants on chromosome 
2 (Chr 2) around Prnp at 128-156 Mbp (Fig. 1 B). 

Of all genes transcribed in pM<I>s, 132 showed nonsynon- 
ymous coding variants between B6.129-Prnp"' r/ "" and B6.129- 
PmpZrM/zrchi ( Fig 1 c all d Table S2). Of these, 62 mapped to 
Chr 2, 4 of which, Sirpa (Matozaki et al, 2009; Chao et al, 
2012), Mertk (Scott et al, 2001), Tyro3 (Seitz et al, 2007), and 
Thbsl (Gao et al., 1996), had been previously described to 
impact phagocytosis (Fig. 1 C). 

Of these four candidate genes, Sirpa (Fig. 1 C) encodes 
the transmembrane SIRPa, is closely linked to Prnp (-2.2 
Mbp), is highly polymorphic in laboratory mice (Sano et al., 
1999), and, uniquely among the four candidates, displays poly- 
morphisms with robust modulatory effects on phagocytosis 
(Takenaka et al., 2007; Legrand et al., 2011; Strowig et al., 2011; 



Theocharides et al., 2012;Yamauchi et al., 2013). Binding of 
CD47 on cognate cells to SIRPa on macrophages transmits a 
"don't eat me" signal that protects the targets from phagocy- 
tosis (Matozaki et al., 2009; Chao et al., 2012). Down-regulation 
of CD47 cell surface expression on ageing or apoptotic cells en- 
hances their phagocytosis (Barclay and Brown, 2006; Matozaki 
et al., 2009). Polymorphisms of the Ig-like amino-terniinal 
domain of SIRPa can affect CD47 binding and modulate 
phagocytosis (Takenaka et al., 2007; Legrand et al., 2011; 
Strowig et al., 2011;Theocharides et al., 2012). These charac- 
teristics led us to prioritize Sirpa as the most plausible candi- 
date for modifying the phagocytic phenotype described in 
B6A29-Pmp z " llI/z " h ' versus Bb.\29-Prnp l "" M mice. 

Sirpa 129 segregates with Prnp~ in seven Prnp~/~ lines 

We then analyzed Sirpa alleles in PfnpZrchi/Zrchi anc [ m s ^ 
additional, independently generated mouse lines carrying tar- 
geted Prnp deletions (Table 1, Fig. 2, and Fig. SI). Restriction 
fragment length polymorphism (RFLP) analysis and DNA 
sequencing revealed that the targeted Prnp allele co-segregated 
with the Sirpa allele of the 129 strain (Sirpa 129 ) in PmpZrchi/Znhi 
mice that had been backcrossed >17 generations to BALB/c 
(C) or for >12 generations to B6 mice in our colony (Fig. 2). 
We also observed co-segregation of Prnp and Sirpa alleles in 
Prnp Zlll ' I/ZrchI mice backcrossed for 12 generations to B6 and in 
p rn pEtti>g/E,ib x backcrossed for <18 generations to B6 mice from 
four different laboratories (Figs. SI and S2 and not depicted). 
Similar results were observed for alleles at Mertk, Tyro3, and 
Thbsl, with the exception of an independently generated 
congenic Pmp~'~ line called "Ngsk" (Table 1 and Table S3). 
Congenic B6.129-Prnp Nxsk/Ngsk mice have the genotypes 
Sirpa 129 ' 129 Mertk 129 " 29 Tyro3 B6/B6 Thbsl B6,B6 , likely reflecting 
a recombination event during repeated backcrossing to B6 
mice. Thus, Sirpa and/or another of the three genes could af- 
fect the hyperphagocytosis phenotype reported in p m pZrdii/z r chi 
mice (de Almeida et al., 2005). 

Functional analysis excludes Prnp and points to Sirpa 
as the gene controlling phagocytosis 

We then sought to determine the contribution of polymor- 
phic Prnp-flanking genes to the phagocytic phenotype of 
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Figure 1. RNA sequencing identifies 
macrophage-expressed genes potentially 
controlling phagocytosis in Prnp Zrchl/z " :hl 
versus Prnp wt/wt mice. (A) mRNA heat map 
highlighting 305 differentially expressed 
genes between B6.129- p r np Zrch " Zr ' :hl and 
B6.129-Pmp»"/ | * t (pM*s, three mice per 
group), ordered by chromosomal location 
(vertical) and unsupervised hierarchical 
clustering (horizontal). Chr 2 (in yellow) 
contains significantly more differentially 
expressed genes (P = 0.001 1, based on 
Z score > 44 and corrected for chromosome 
size). Individual genes are listed in Table S1. 
(B) Physical distribution of SNPs unique to 
the B6.129-Pmp Zrc W ZrcW pM4> transcriptome 
clustered on Chr 2 (in yellow; two-tailed 
Student's ttest, P = 0.019) at 128-156 
Mbp. (C) Filtering strategy to prioritize 
genes controlling phagocytosis in Prnp Zrchl/ 
Zrchl versus Pmp wt,wt mice. Sequential filters 
were applied to a total of 22,981 SNPs dif- 
fering between the B6A2S- Prnp Zrch " Zr < :hl and 
B6.129-Pmp ,vt '' wt pM$ transcriptomes. SNPs 
affecting protein coding sequence included 
nonsynonymous SNPs and insertions/dele- 
tions. The asterisk indicates phenotypic 
definition based on gene ontology annota- 
tions and PubMed searches. Sirpa, Mertk, 
Thbsl, and Tyro3 display nonsynonymous 
SNPs between B6^29-Prnp z ' ch " z ' c> > 1 and 
B6.129-Pmp"^ 1, "mice, reside on Chr 2, and 
are involved in phagocytosis. However, only 
Sirpa (purple) is known to display polymor- 
phisms that modulate phagocytosis. 
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Pmp" 7- mice.The phagocytic activity of primary macrophages 
was compared between Prnp wt/u " and three different strains of 
Prnp~ /- mice differing at all four, two, or none of the closely 
linked polymorphic genes. When offered apoptotic thymo- 
cytes as targets, B6A29-Pmp ZrchI/ZrchI pM$s displayed higher 
rates of phagocytosis than those from heterozygous B6.129- 
p mp wt/z,chi orWT B6A29-Prnp u " /u " congenic littermates (Fig. 3, 
A and B, left) in agreement with a previous study (de Almeida 
et al., 2005). We also observed this hyperphagocytic phenotype 
m BMDMs of B6A29-Pmp ZrchI/z,M and B6A29-Prnp N s sk/N v k 
mice (Fig. 3 B, middle). In contrast, pM$s from \29-Prnp m z' Edb s 



mice (Fig. 3 A, right; Manson et al., 1994) displayed similar 
phagocytic activity as 129-Prnp""'"" pM$s (Fig. 3 B, right). 
\29-Prnp E,,t ' x/m ' s were generated from 129 ES cells and main- 
tained as a co-isogenic 129 line (Fig. 3 A, right). The two 
strains showed no differences in microsatellite markers linked 
to Pmp (Fig. 3 C). Consequently, ablation of Prnp was insuffi- 
cient to drive hyperphagocytosis, suggesting the influence of 
additional genetic elements in Zrchl and Ngsk mice. 

To define the 129-derived regions flanking Pmp, we per- 
formed microsatellite analyses of Chr 2 in B6A29-Prnp~ / ~ 
mice of the Zrchl and Ngsk lines. B6A29-Pmp ZrM/ZrM mice 
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Table 2. Summary statistics of RNA sequencing 



Mouse 


Total reads Mapped reads 


Mean coverage 


# SNPs 


# Het. SNPs 


# Horn. SNPs 


# Indels 


# Nonsyn. 
variants 


# Nonsyn. 
variant genes 


WT 1 


58,838,924 56,801,099(96.5%) 


58.6x 


102,465 


1 4,249 


88,216 


6,860 


2,902 


362 


WT2 


56,971,575 55,447,441 (97.3%) 


60. 1x 


102,756 


12,591 


90,165 


7,451 


2,874 


361 


WT 3 


77,632,102 75,325,378 (97.0%) 


67.4x 


117,354 


15,501 


101,853 


8,356 


3,360 


415 


All WT 






84,205 


9,535 


74,670 


5,644 


2,277 


197 


KO 1 


71,507,393 69,906,771 (97.8%) 


64.6x 


117,979 


15,146 


102,833 


8,164 


3,517 


466 


K0 2 


70,456,814 68,728,168 (97.5%) 


62.7x 


114,906 


14,943 


99,963 


8,185 


3,391 


466 


K0 3 


75,295,650 73,039,778 (97.0%) 


65.4x 


117,592 


15,539 


102,053 


8,405 


3,515 


488 


All KO 






92,176 


9,569 


82,607 


6,397 


2,682 


264 



WT, B6.129-Pmp""/""; K0, B6.129- p r np Zrch V z ' chl ; #, number; Het, heterozygous; Horn., homozygous; Indels, insertions/deletions; Nonsyn., nonsynonymous. 



harbored 129-derived genomic material in a region of >48 
Mbp flanking Prnp on Chr 2 (including Sirpa, Mertk, Tyro3, 
and Thbsl alleles) despite >12 generations of backcross to B6 
(Fig. 3 C and Fig. S2). 

Compared with B6A29-Prnp Zn ' l[/ZrM , B6.129-Pmp N ^ /N ^ 
mice contained a smaller 129-derived interval on Chr 2 
(Fig. 3 C). These mice were homozygous for 129 alleles at 
Sirpa and Mertk and homozygous for B6 alleles at Tyro3 and 
Thbsl (Fig. 3 C). Because B6.129-Pmp N ^ /f ** mice retained 
the hyperphagocytic phenotype, variations in Tyro3 or Thbs 1 
alleles did not explain the differential phagocytic activity 
of Prnp u " /u " and Prnp"'" mice, excluding these two genes as 
candidates (Fig. 3 C). 



We identified only one nonsynonymous single nucleotide 
polymorphism (SNP) in Mertk 129 versus Mertk M alleles (Chr2: 
128776396, Att/Gtt). This SNP, previously reported in CD-I 
mice (SNP accession no. rs27446326 from dbSNP Build 137), 
results in a conservative isoleucine-to- valine substitution within 
the transmembrane region (15 16V) and is predicted by the 
SIFT algorithm (Ng and Henikoff, 2001; Kumar et al., 2009) 
to be tolerated and unlikely to impact Mertk function. All dis- 
ease-causing MERTK mutations identified so far in patients 
suffering from retinitis pigmentosa (Gal et al., 2000), Leber 
congenital amaurosis (Li et al., 2011), or cancer (Greenman 
et al., 2007) were located in the extracellular or cytoplasmic 
but not in the transmembrane region. Moreover, one natural 
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Figure 2. Sirpa allelotype of the Prnp' 1 ' lines 
analyzed in this study. (A) RFLP analysis of Sirpa 
allelotypes. The Sirpa' 29 allele segregates with the 
Prnp" allele in Prnp" 1 " mouse lines stemming from 
seven independent targeting events (see Table 1). The 
8th through 1 1th lanes show respective reference 
strain DNA. (B) The red box indicates protein se- 
quence alignment of mouse SIRPa Ig-like variable 
domain (IgV) confirming the 129 Sirpa allelotype in 
each one of the seven Prnp" 1 " lines shown in A. The 
blue box indicates SIRPa protein reference sequences 
for the B6, 129/Sv, and 129/Ola strains. Polymor- 
phisms are highlighted in yellow. For each group one 
mouse was analyzed. 
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Figure 3. Ablation of Prnp is not sufficient to enhance phagocytosis of apoptotic cells. (A) Generation of B6.129 congenic (left) and 129 
co-isogenic (right) Prnp~/~ lines. In B6.129 congenic mice, residual 129 genomic material derived from the ES cell (brown) is inevitably present, 
particularly in the region flanking the targeted locus. (B, left) pM<t>s from congenic B6.M9-Prnp Zrchi/Zrchl mice showed higher phagocytic activity 
than those from B6.1 29-Prnp wt/w < (set as 100%) and B6.1 29- Prnp M/Zrchl littermates. Blue, black, and gray (here and henceforth) indicate data from 
independent experiments, (middle) B6.1 29-Prnp Zrch " Zrch! and B6.129-Prnp N ^ N s sk BMDMs showed higher phagocytic activity than Be-Prnp"""' mac- 
rophages (set as 100°/o). (right) The phagocytic activity of 1 29-Prnp Edb9/Edl ' s and 1 29-Prnp wt/r " macrophages was similar, even when two different 
macrophage/thymocyte (M:T) ratios were used (1:10 and 1:1; one-way AN OVA, Bonferroni's multiple comparisons test). Mean 1 29-Prnp w,/wt 
phagocytic rate at 1:10 M:T ratio was set as 100%. n.s, not significant; **, P < 0.005; ***, P < 0.001 ; ****, P < 0.0001. Error bars indicate SD. (C) STR 
analysis of Chr 2. Name and physical position of each STR marker (colored boxes) are indicated on the left. B6.1 29-Prnp ZrM/z ' M and B6.129- 
p rn pNgsk/Ngski m | ce contained 129-derived genetic material flanking Prnp, whereas the markers flanking Prnp in 129-Pmp M, f /& *!' and 129- Prnp wt/wt 
mice were identical. Data show representative mice of at least three animals analyzed per group. 
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Figure 4. Sirpa 86 allele segregates with enhanced inhibition of phagocytosis in the absence of Prnp. (A) Breeding scheme applied to generate 
B61 29-Prnp Zrchl ' Zrchl mice containing Sirpa 86 and Sirpa' 29 alleles. B61 29-Prnp Zrch " Zrchl Sirpa B6 ' B6 tga2ffs/'g m j C e (F 0 ) were backcrossed to B61 29-Pmp z ' M l z ' M 
Sirpa 129 ' 129 mice for two generations to produce B6M9-Prnp Zrchl/Zrchl littermates, with all combinations of Sirpa m and Sirpa' 29 alleles (F 3 ) used to assess 
phagocytic activity as in C. (B) RFLP analysis of Sirpa allelotypes (B6 vs. 129). Surprisingly, fga20 mice (first through third lanes) displayed combinations 
of Sirpa 86 and Sirpa 129 alleles. The fourth through sixth lanes show B^29-Prnp Zrchl ' Zrchl mice labeled in A (red arrows). Controls are reference DNA. All 
animals entering the study were analyzed. (C) Phagocytic hyperactivity of B612<3-Prnp Zrchl/z " :hl pM<I>s was associated with the Sirpa 129 " 29 allelotype. Data 
from three independent experiments (blue, black, and gray) normalized against mean BG129-Prnp Zrch " Zrchl Sirpa B6/Be phagocytic rates. One-way AN0VA, 
Bonferroni's multiple comparisons test: * P < 0.05; ***, P < 0.001. Error bars indicate SD. (D) STR analysis documenting the boundaries of B6 versus 129- 
derived genome in Chr 2 of B6129-Pmp Zrc '' ,/Zrc '''with different combinations of Sirpa alleles (B6 vs. 129). Name and position of each STR marker (colored 
box) on Chr 2 are indicated on the left. Data show representative mice of at least three animals analyzed per group. 



non— disease-causing variant in the transmembrane region has 
been reported in humans (15 18V; Gal et al., 2000; Li et al., 
2011), which corresponds to the I516V variant in our mice, 
again suggesting that the variant identified in our mice is un- 
likely to impact Mertk function. Collectively, these findings 
suggested that Sirpa 129 homozygosity conferred differential levels 
of phagocytosis in otherwise genetically homogeneous mice and 
that the presence of at least one Sirpa 86 allele reduced rates of 
phagocytosis compared with Sirpa 129 homozygotes. 

Sirpa 86 association with inhibition of phagocytosis 

We observed that some of our PrP c -overexpressing tgz20 
transgenic mice on a mixed B(y\29-Prnp ZrdlI ' z " :hI background 
(Fischer et al., 1996) were segregating Sirpa 86 and Sirpa 129 
alleles, likely reflecting an unintentional cross with B6-Prnp""'"" 
breeders in the early days of this colony (unpublished data) .We 
crossed m\29-Prnp ZrchI,ZrM Sirpa M/B6 tgz.20^ mice to B6129- 
PrnpZrAi/zM Sirpa um29 mlce pro d uc ing B6129-Prnp ZrM/z,M 
littermates with all three combinations of Sirpa 86 and Sirpa 129 



alleles (Fig. 4, A and B). pM$s from B6l29-Prnp Zrcl, ' /ZlThI 
Sirpa 129 ' 129 mice showed higher rates of phagocytosis (Fig. 4, 
C and D). Hence, compared with Sirpa 1291129 , a single Sirpa 86 
allele reduced phagocytosis rates in the absence of Prnp. These 
results supported the hypothesis that homozygosity for Sirpa 129 
rather than Prnp deficiency was responsible for the hyperpha- 
gocytic phenotype of Prnp Zrd ' I/ZrchI mice. 

Inhibition of phagocytosis segregated 
with Sirpa 86 in Prnp wt/wt mice 

Sirpa alleles and Prnp ablation may impact phagocytosis in- 
dependently of each other. We therefore asked whether dif- 
ferent Sirpa alleles could modify phagocytosis in Prnp""'"" 
mice.We genotyped 103 pups from B6.l29-Pmp"'' /ZnhI Sirpa 86 ' 129 
intercross mice for recombination between Prnp and Sirpa and 
observed two events producing a Prnp"' 1 -Sirpa 129 and a Prnp z " hl - 
Sirpa 86 haplotype (Fig. 5 A). Mice carrying the original and the 
recombinant haplotypes were then bred to generate congenic 
mice with different Prnp-Sirpa haplotypes (Fig. 5,B and C). 
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Figure 5. Hyperphagocytosis is associated with the Sirpa U9/129 allelotype but not with Prnp gene dosage. (A) Breeding scheme to generate 
recombinant congenic B6.129 mice with different combinations of Prnp (WT and Zrchl) and Sirpa (B6 and 129) alleles. BS.12<3-Prnp r " /Zrch! mice (F,) were 
intercrossed, and F 2 offspring were screened for meiotic recombination between Prnp and Sirpa. 103 F 2 mice with different combinations of Prnp geno- 
types and Sirpa allelotypes were obtained. Two mice (red frames) were found to carry the recombinant haplotype, Prnp w '-Sirpa' 29 and Prnp z ' chL Sirpa B6 
and were interbred to generate B5A29- Prnp Zrchl/Zrchl and B6.129-Pmp"" / "" mice homozygous for Sirpa 86 or Sirpa' 29 alleles. (B) RFLP analysis to discriminate 
between Sirpa 129 and Sirpa 86 . BGM9-Prnp wt/w > (first, third, and fifth lanes) and BGA29-Prnp Zrchl/Zrchl (second, fourth, and sixth lanes) mice with different 
combinations of Sirpa 129 and Sirpa Be alleles. Controls are reference DNA. All animals entering the study were analyzed. (C) Protein sequence variants (yel- 
low) in SIRPa Ig-like variable domain (IgV). Homozygous recombinant congenic B6.129 mice displayed four combinations of Prnp-Sirpa haplotypes (black 
boxes). Blue boxes show SIRPa sequence of reference B6 and 129/Sv strains. For each group one mouse was analyzed. (D) Sirpa' 29/ ' 29 , but not Prnp^^, 
was associated with hyperphagocytosis in pM<t>s from congenic B6.129 mice. Data are from two independent experiments (blue and black); phagocytosis 
rates were normalized against B6.129 Prnp wt/w < Sirpa B6/B6 . One-way AN0VA, Bonferroni's multiple comparisons post-test:*, P < 0.05; ****, P < 0.0001. Error 
bars indicate SD. (E) Recombination breakpoints between Prnp and Sirpa in congenic B6.129 mice. For each STR marker (colored box), name and position 
on Chr 2 are indicated on the left. Data show representative mice of at least three animals analyzed per group. 



The various Prnp-Sirpa haplotypes did not affect the fre- 
quency of macrophages in peritoneal lavages, nor did Sirpa 
allelotypes affect PrP c expression (Fig. 6, A— C). Also, we did 



not observe differences in Sirpa rnRNA levels between pM<I>s 
from m.\29-Pmp Zrd ' I,z " hI and B6.129-Pmp"" / "" mice (375.60 
vs. 356.34 fragments per kilobase of transcript per million 
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mapped reads [FPKM], respectively, P = 0.73) or SIRPct pro- 
tein levels in mice with different Prnp-Sirpa haplotypes (Fig. 6 D). 
These results excluded any influence of Prnp or Sirpa alleles on 
the expression of each other. 

In congenic B6.129 mice with different combinations of 
Prnp and Sirpa alleles (WT and ZrchI,B6 and 129, respectively), 
hyperphagocytosis segregated with Sirpa 129 ' 129 irrespectively of 
the Prnp status (Fig. 5, D and E). We compared phagocytosis by 
BMDMs isolated from congenic C.\29-Prnp z " l '" z ' d '' mice with 
different combinations of Prnp and Sirpa alleles (Fig. 7, A— C). 
In line with our findings with congenic B6.129 mice, a higher 
phagocytosis rate correlated with Sirpa 129 ' 129 independently of 
the Prnp alleles (Fig. 7, D and E). Together, these results exclude 
a role for PrP c and support Sirpa polymorphisms as a major 
regulator of phagocytosis in mouse strains carrying either WT 
or deleted Prnp alleles. 

DISCUSSION 

In the 1990s the Prnp gene became an extremely popular 
target among mouse geneticists and rapidly advanced to one 
of the most frequently inactivated elements of the mouse 
genome. As a result, we have indeed identified at least seven 
published mouse lines derived from independent gene- 
targeting events, all of which were analyzed in the present 
study. Remarkably, we discovered that all of these lines suffer 
from the same set of confounders when kept as congenic: 
a 129-derived genomic region sporting functionally impact- 
ful polymorphisms. 

The alleged control of phagocytosis of apoptotic cells by 
Prnp was critically questioned when we found that, in con- 
trast to congenic B6.129- or C.\29-Pmp z " wz " hI and B6.129- 
p rn pNg S k/N gs k micC; co _i SO genic 129-Prnp Edb s /E < lh £ mice did not 
display a hyperphagocytic phenotype. Moreover, microsatel- 
lite-assisted breeding of B6129-Prnp z mice revealed 

the existence of a 129-derived Chr 2 segment which, when 
homozygous, correlated with a higher rate of phagocytosis 
independent of Prnp. These data suggested that Prnp does not 
impact phagocytosis and that such conclusions in previous 
studies resulted from unrecognized polymorphic flanking 
genes (de Almeida et al., 2005). Detailed genomic, transcrip- 
tomic, and functional analysis of candidate genes in the intro- 
gressed 129-derived region in 'B6.129-Prnp z " :l ' I ' ZrchI mice 
supported a focus on the well-characterized inhibitor of 
phagocytosis Sirpa, whose polymorphisms alter this pheno- 
type (Takenaka et al., 2007; Legrand et al., 2011; Strowig 
et al., 2011). Indeed, Sirpa 129 , rather than Prnp~, segregated 
with increased phagocytosis. 

B6A29-Prnp Zreh " ZKhI and B6A29-Prnp""' u " macrophages 
displayed sequence variations, or differential expression, of 
several other genes in addition to Sirpa. In particular, Mertk 
(Scott et al., 2001), Tyro3 (Seitz et al, 2007), and Thbsl (Gao 
et al., 1996), which participate in phagocytosis and are linked 
to Prnp, were polymorphic between the two strains. However, 
the genetic and functional considerations detailed above sug- 
gested that they were unlikely candidates. 
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In our study, we have used mouse genetics to investigate 
the role of PrP c and its polymorphic flanking genes in modu- 
lating phagocytosis of apoptotic cells. Compared with other 
experimental approaches, including gene silencing and over- 
expression of exogenously introduced genes, formal genetics 
has the advantage of assessing the differential impact of physi- 
ological versus absent expression of the gene under scrutiny. 
Moreover, silencing or overexpressing Sirpa, or other modu- 
lators of phagocytosis, could alter the delicate balance between 
pro- and antiphagocytic signals (Hochreiter-Hufford and 
Ravichandran, 2013) and impact phagocytosis independently 
of whether this modulator is responsible for the difference 
between Prnp~'~ and Prnp +n macrophages. 

Enzymatic removal of PrP c from macrophage surface 
through phosphoinositol phospholipase C was reported to re- 
sult in an increased phagocytosis of apoptotic cells, seemingly 
confirming the role of PrP c in this process (de Almeida et al., 
2005). However, this treatment genetically affects all glyco- 
sylphosphatidylinositol-anchored proteins, many of which are 
involved in the control of phagocytosis (Dorahy et al., 1996; 
Cebecauer et al, 1998; Oliferenko et al., 1999; Kellner-Weibel 
et al., 2000; Peyron et al., 2000; Chimini, 2001; de Almeida and 
Linden, 2005). Therefore, it is questionable whether this effect 
is directly mediated by loss of surface PrP c . 

Our study focused on the presumed role of PrP c in 
phagocytosis of apoptotic cells. Phagocytosis of zymosan par- 
ticles and latex beads were also reported to be influenced by 
Prnp using primary macrophages or immortalized cell lines 
derived from Prnp~'~ and Prnp +n mice (de Almeida et al., 
2005; Nitta et al, 2009; Uraki et al, 2010). Our results raise 
the question of whether Sirpa polymorphisms may affect the 
responses to zymosan or phagocytosis of latex beads as well. 

The co-segregation of Sirpa 129 with Prnp~ in all the 
Prnp~'~ hnes examined has implications beyond the control 
of phagocytosis. Because the Sirpa sequence of most mouse 
strains used for backcrossing of Prnp~'~ mice (B6, BALB/c, 
and FVB) differs from 129 mice (Sano et al., 1999; Takenaka 
et al., 2007), Sirpa is a potential confounder in most studies of 
Prnp - ' - mice. Importantly, the expression pattern of Sirpa 
mirrors that of Prnp and includes myeloid cells, neurons, and 
pancreatic (3 cells (Linden et al, 2008; Matozaki et al, 2009). 
Finally, PrP c -deficient mice share numerous subtle pheno- 
types with S/Vpa-ablated mice, including neuroimmunologi- 
cal (Tomizawa et al., 2007; Tsutsui et al., 2008), behavioral 
(Nico et al, 2005; Ohmshi et al, 2010; Gadotti et al, 2012), 
and metabolic phenotypes (Kobayashi et al., 2008; Strom 
et al., 2011). These observations suggest that additional pheno- 
types attributed to Prnp may actually reflect the effects of 
Sirpa polymorphisms. 

Based on studies with Pmp~'~ mice, PrP c has been re- 
cently reported to either induce or inhibit autophagy (Oh 
et al, 2008, 2012; Korom et al, 2013) and to be indispensable 
or dispensable to mediate A(3 toxicity (Lauren et al., 2009; 
Balducci et al., 2010; Calella et al, 2010; Gimbel et al, 2010; 
Kessels et al., 2010; Cisse et al., 2011; Larson et al., 2012; Um 
et al., 2012). At least some of these controversies could reflect 
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Figure 6. Sirpa allelotype does not influence pM<I> abundance or PrP c and SIRPa expression. (A) Peritoneal lavages were gated for nucleated 
cells (left), and the percentage of F4/80+ cells was assessed (middle panel, right peak), (right) Overlay of histograms of F4/80 (dark blue), isotype antibody 
(light blue), and unstained peritoneal cells (gray) confirming the specificity of the staining. (B) Percentage of F4/80 + macrophages in peritoneal lavage 
cells of experimental mice. The dotted lines delineate individual experiments in which groups of mice (n = 3-10 for each haplotype) were analyzed. Error 
bars indicate SD. (C) PrP c ELISA analysis of pM<I> lysates showed no difference (n.s.) in PrP c levels between B6.129-Pmp vvt/ "' t mice with different Sirpa 
allelotypes. Two-tailed unpaired Student's f test: P = 0.241 ; n = 3-4 for each haplotype. n.d., not detectable. (D) Unchanged expression of the SIRPa 
cytoplasmic tail (CT) in cultured pM<l> lysates as quantified to actin levels. SHPS-1 K0 macrophage lysates lacking the cytosolic tail of SIRPa were used 
as negative control. One-way AN0VA: n.s., not significant; P = 0.27; n = 3-4 for each haplotype. Horizontal bars indicate mean. 



the effect of polymorphic Prnp-flanking genes, and these 
findings merit critical reconsideration in light of our data. 

The risk of flanking-gene problems can be reduced by 
reversing the phenotype of gene-targeted mice through trans- 
genetic restoration of gene expression (Wolfer et al., 2002). 
With Prnp, this was often done by using PrP c -overexpressing 
tga.20 mice, which were generated from PffipZ^i/znhi lrLlce anc J 
subsequently distributed to many laboratories. However, we 

JEM Vol. 210, No. 12 



found that the tgz20 colony was contaminated with multiple 
combinations of Sirpa alleles early after its establishment. 
Hence the rescue of Prnp deficiency through a tga.20 allele 
does not necessarily prove that the observed changes in pheno- 
type are PrP c dependent. 

It has been long suspected that genes flanking targeted 
loci may confound the interpretation of KO experiments 
(Smithies and Maeda, 1995; Gerlai, 1996). Recommendations 
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Figure 7. Phagocytic activity level is associated with Sirpa 1291129 allelotype but not with Prnp gene dosage also in C.129 congenic 
mice. (A) Breeding scheme to generate recombinant congenic C.129 mice with different combinations of Prnp (WT and Zrchl) and Sirpa (C and 
129) alleles. C\29-Prnp r " /Zrch ' mice (F,) were intercrossed, and the occurrence of meiotic recombination between Prnp and Sirpa was assessed in 
the F 2 offspring. 243 F 2 mice with different combinations of Prnp genotypes and Sirpa allelotypes were obtained. Four mice (red frames) were 
found to carry a recombinant haplotype, and selected ones were bred to generate CA29-Prnp Zrchl/Zrchl and C.129-Prnp l,rt/ "' t mice homozygous for 
Sirpa 0 or Sirpa' 29 alleles. (B) RFLP analysis used to discriminate Sirpa' 29 from Sirpa c alleles. C.1 29-Prnp wt/wt (first, third, and fifth lanes) and 
CM°)-Prnp z ' ch ' lz ' M (second, fourth, and sixth lanes) mice with different combinations of Sirpa 129 and Sirpa c alleles. Controls are reference DNA. All 
animals entering the study were analyzed. (C) Protein sequence alignment of mouse SIRPot Ig-like variable domain (IgV) illustrates recombinant 
congenic C.129 mice with different combinations of Prnp and Sirpa' 29 versus Sirpa c (gray border). Blue boxes are SIR Pa protein reference se- 
quences for the B6, 129/Sv, and C strains (yellow: polymorphisms). For each group one mouse was analyzed. (D) Sirpa' 29 !' 29 allelotype, but not the 
absence of Prnp, was associated with hyperphagocytosis of BMDMs in congenic C.129 mice. Data are from three independent experiments (blue, 
black, and gray). Mean phagocytosis rates of C.129 Prnp Zrch '/ Zrchl Sirpa Be/B6 were set as 100°/o. One-way AN0VA, Bonferroni's multiple comparison 
post-test: #, P = 0.075; ****, P < 0.0001. Error bars indicate SD. (E) STR analysis documents the result of trans-allelic meiotic recombination be- 
tween Prnp and Sirpa. Blue dashed line indicates the location of Prnp and Sirpa. For each STR marker (colored box), name and position on Chr 2 
are indicated on the left. Data show representative mice of at least three animals analyzed per group. 
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for breeding schemes and genetic quality control to avoid this 
problem were formulated (Wolfer et al., 2002; Ridgway et al., 
2007; Crusio et al., 2009) but rarely implemented (Holmdahl 
and Malissen, 2012). Consequently, reports of flanking gene 
effects are anecdotal and rely on the demonstration that a 
KO mouse phenotype is lost upon crossing to other strains 
backgrounds and reduction in length of the 129-derived re- 
gion flanking the targeted locus (Kanagawa et al., 2000; de 
Ledesma et al., 2006;Eisener-Dorman et al., 2010). Our find- 
ings provide an alternative mechanism underlying prior spe- 
cific claims in the prion field and strengthen the requirement 
for rigorous genetic, genomic, and functional analyses to sup- 
port the assignment of specific phenotypes to targeted muta- 
tions in mice. 

MATERIALS AND METHODS 

Mice. The following mice were analyzed: Prnp Zn ' l ' l/Zn ' 1 ' 1 mice (Biieler et al., 1992) 
on a mixed B6129 background or backcrossed to B6 for >12 generations or to 
C for > 17 generations and congenic WT B6 and C mice purchased from Harlan 
or Charles River or bred in house; Pmp Nssk/Nssi! mice (Sakaguchi et al., 1995) 
extensively backcrossed to B6; Pmp Edh £ /Edh& mice on a pure 129/Ola back- 
ground (Manson et al., 1994) and co-isogenic WT 129/Ola mice; PmpZr&n/zxhn 
mice (Rossi et al., 2001) on a mixed B6129 background; Prnp GFF/GFP mice 
(Heikenwalder et al., 2008) backcrossed to B6 for 10 generations; tga20 mice 
(Fischer et al., 1996) on a mixed B6129 p m f- rM/ZrM background; and SHPS-1 
(Inagaki et al., 2000) backcrossed to B6 for >5 generations. Genotype of mice 
was verified as indicated in the original description of each line (Table 1). To 
introduce Sirpa B6 into B6129-Pmp ZrcllI/ZrchI mice, B6\29-Pmp ZrcWZrchl Sirpa B6/B6 
tgz2Qs/'s nuce were crossed with R6\29-Pmp ZrcWZrchI Sirpa 129/ ' 29 mice (Fig. 4 A). 
To generate congenic 'B6.V29-Pmp ZrihI/ZKl ' 1 and Pmp u ' t/li,i mice with all possible 
combinations of Sirpa 36 and Sirpa 129 alleles, B6.129-Pmp"'' /Zr,-W Sirpa B6/l29 mice 
were intercrossed, and recombination events between Prnp and Sirpa loci were 
identified by Prnp genotyping and RFLP analysis of Sirpa exon 2 in the offspring 
(Fig. 5, A and B). Mice with a recombinant haplotype (referred to as obligate 
carriers) were further crossed. The same approach was followed to generate 
congenic Q.\29-Pmp Zfcb!/Zfd,! and Prnp wt/m mice with all possible combina- 
tions of Sirpa c and Sirpa 129 alleles (Fig. 7, A andB). All animal experiments were 
performed in compliance with the Swiss Animal Protection Law, under the 
approval of the Veterinary office of the Canton Zurich. 

Allele discrimination analysis. To distinguish strain-specific alleles, 
PCR-RFLP analysis (Sirpa, Tyro3, and Thbsl) or sequencing (Mertk) were 
performed using primers and restriction enzymes as detailed in Table S4. 

Phagocytosis assay with pM$s. On day 0, peritoneal cells were collected 
from age- and gender-matched adult mice according to standard methods 
(Fortier and Falk, 2001), resuspended to 2.2 x 10 6 in 5 ml DMEM (Gibco) 
with 10% FBS (Gibco), 1% GlutaMAX (Gibco), and antibiotics, plated in 
one-well slide chambers (Falcon; BD), and incubated at 37°C. After 2 h, 
nonadhering cells were removed and remaining cells (referred to as cultured 
pM<I>s) were kept overnight in culture. On day 1, thymocytes were harvested 
from 3— 8-wk-old mice of the same background as the genotypes under in- 
vestigation, incubated for 3 h at 37°C in the presence of 1 uM dexamethasone 
to induce apoptosis, washed, and added to the cultured pM^s at different ra- 
tios for 1 h. For each experimental group, technical replicates (at least dupli- 
cates) were performed. Cells were washed, fixed with formalin, and stained 
with May-Griinwald-Giemsa. Slides were randomized and scanned with 
NanoZoomer (Hamamatsu Photonics), random fields (50 fields/slide) were 
generated, and phagocytosis rate was assessed using in-house— developed soft- 
ware. The same operator blinded to experimental groups validated all soft- 
ware-analyzed slides. For each technical replicate, >800 cells were evaluated 
(indicated as a dot in the graphs). Phagocytosis rate was defined as the per- 
centage of macrophages having at least one apoptotic body. 



Phagocytosis assay with BMDMs. On day 0, femurs were collected from 
age- and gender-matched adult mice and flushed with RPMI1640 medium 
(Gibco) containing 10% FBS (Gibco), 1% GlutaMAX (Gibco), and antibiot- 
ics. BM cells were plated into 6-well plates at 2 X 10 6 /well in 3 ml of culture 
medium containing 10 ng/ml of macrophage colony stimulation factor (In- 
vitrogen) and cultured overnight at 37°C. On day 1, cells were transferred to 
new 6-well plates to remove loosely adherent stroma-like cells. On day 6, 
cells were harvested, adjusted to 5 x 10 5 cells in 500 ul, and plated into 24- 
well plates. Apoptotic thymocytes (see above) were suspended at 10 6 /ml in 
PBS and labeled with 20 ng/ml of the pH-sensitive dye pHrodo red, SE 
(Invitrogen) for 30 min at room temperature and washed, and 500 ul of cell 
suspension was added to BMDM culture for 1 h at 37°C. After washing, 
BMDMs were harvested with Accutase (Invitrogen) and gentle scraping and 
stained with FITC-labeled anti-CDllb or isotype control antibody (BD), 
and rate of phagocytosis was analyzed by flow cytometry as the percentage of 
pHrodo positivity among CDl lb + cells. At least 10,000 events were acquired 
in the living gate. 

Genome-wide and Chr 2 short tandem repeat (STR) analysis. 

Genome-wide STR analysis was performed as previously described (Bremer 
et al., 2010). In brief, purified genomic DNA (gDNA) was amplified using 
multiplex PCR with nuorescently labeled primers (FAM, NED, and VIC; 
Applied Biosystems), diluted, denatured, and subjected to capillary electro- 
phoresis on a 3130.\'/ Genetic Analyzer (Applied Biosystems). Various mouse 
strains were used for calibration. Peak detection, binning, and allele calling 
were performed using the GeneMapper software (Applied Biosystems) in 
combination with in-house— developed software. 

High-resolution STR genotyping on Chr 2. Novel microsatellite mark- 
ers were generated as previously described (Ivakine et al., 2006; primer se- 
quences in Table S4). DNA was amplified for 35 cycles: 30 s at 94°C, 30 s at 
55°C, and 1 min at 72°C using Multiplex Master Mix (QIAGEN).The PCR 
product was diluted 1/20, and 1 \jX mixed with a 10-|-il mixture containing 
formamide and a 500LIZ size standard. Samples were resolved on the 3730x/ 
DNA Analyzer (Applied Biosystems). Alleles were sized, in comparison with 
standards, by viewing electropherograms in GeneMapper. 

Library preparation for RNA sequencing. Total RNA was isolated from 
cultured pM^s using the RNeasy mini kit (QIAGEN), snap frozen, and kept at 
— 80°C until further analysis. RNA quality was determined with Qubit (1.0) 
Fluorometer (Life Technologies) and Bioanalyzer 2100 (Agilent Technologies). 
Only samples with a 260 nm/280 nm ratio between 1.8 and 2.1 and a 28S/18S 
ratio between 1 .5 and 2 were further processed. The TruSeq RNA Sample Prep 
kit v2 (Illumina) was used in the succeeding steps. In brief, total RNA samples 
(1 ug) were poly A enriched and reverse transcribed into double-stranded 
cDNA. TruSeq adapters were ligated to double-stranded cDNA. Fragments 
containing TruSeq adapters on both ends were selectively enriched with PCR. 
Quality and quantity of enriched libraries were validated using Qubit (1.0) 
Fluorometer and Caliper GX LabChip GX (Caliper Life Sciences). The pro- 
duct is a smear with a mean fragment size of °-'260 bp. Libraries were nor- 
malized to 10 nM in Tns-Cl 10 mM, pH 8.5, with 0.1% (vol/vol) Tween 20. 

Cluster generation and RNA sequencing. TruSeq PE Cluster kit 
v3-cBot-HS (Illumina) was used for cluster generation using 2 pM of pooled 
normalized libraries on the cBOT. Sequencing was performed on Illumina 
HiSeq 2000 paired-end at 2 x 101 bp using the TruSeq SBS kit v3-HS (Illumina). 
Sequences have been deposited in NCBI Sequence Read Archive under SRA 
identifier no. SRP029761. 

RNA sequencing data analysis. Reads were quality-checked with 
FastQC. Low-quality ends were clipped (3 bases from the start, 10 bases from 
the end). Trimmed reads were aligned to the reference genome and transcrip- 
tome (FASTA and GIFF files, respectively, downloaded from the UCSC 
mmlO database) withTopHat version 2.0.6.TopHat was run with the follow- 
ing options: mate inner distance set to 30, the corresponding SD to 100, and 
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maximum 10 multi-hits were allowed in the alignment. Reads that did not 
align at the first attempt were split into 25-base-long sections on which a 
second attempt of alignment was performed. In ~B6A29-Pmp z "' hl/Zrchl pM$s, 
all reads mapping to Prnp between BstEII and EcoRI restriction sites in exon 
3 and in the 3' UTR, respectively, possibly representing a fused mRNA con- 
taining neo and residual Prnp, reported in the brain of these mice (Biieler et aL, 
1992), were excluded from the analysis. Polymorphisms were detected using 
GATK version 2.1.8 using the following options: baq Gap open penalty 
(whole-genome analysis) set to 30; minimum consensus coverage to geno- 
type indels set to 8 (default: 5); minimum base quality score and minimum 
variants phred score set to 15. Variants were annotated using snpEFF version 
3.0, and distribution of the reads across genomic isoform expression was 
quantified using RSEM (Li and Dewey, 2011). Cufflinks version 2.0.2 and 
differentially expressed genes listed with its utility Cuffdiff using default op- 
tions were used. 

All remaining data and statistical analyses, formatting, and picture gener- 
ating were produced via in-house R-scripts (R version 2. 15. 2). The SIFT algo- 
rithm with default settings was used to predict the impact of a nonsynonymous 
SNP in Mertk (Kumar et al., 2009). 

Strpa open reading frame sequencing. Cultured pM<I> cDNA and 
gDNA were used to PCR amplify and sequence coding and untranslated re- 
gions of mouse Sirpa as previously described (Takenaka et al., 2007). Primers 
(Micro synth) are listed in Table S4. For Pmp Rkn/R3m mice, sequencing was per- 
formed on gDNA obtained from Institute of Physical and Chemical Re- 
search, Japan, and from cDNA isolated from Pmp Rhl/Iik "— derived hippocampal 
like cells (HpL3-4; Kuwahara et al., 1999). For p mp ^/Rcm0 miCG) only gDNA 
from one single mouse was available, and sequencing was performed only on 
the highly polymorphic exon 2 of Sirpa. 

"Western blotting. Western blotting on pM<I> cell lysates was performed ac- 
cording to standard methods using polyclonal anti-S/'r/M cytoplasmic tail 
(CT) antibody (IMGENEX). 

PrP c -ELISA. PrP c was quantified in cultured pM4> lysates by sandwich ELISA 
as described previously using POM1 and POM2 antibodies {Polymenidou 
et al.,2008). 

Flow cytometry analysis. Freshly isolated peritoneal cells were incubated 
with FITC-labeled anti-F4/80 or isotype control antibody (BD). Thymocyte 
apoptosis was assessed using the FITC AnnexinV Apoptosis Detection kit II 
(BD). Flow cytometry was performed using a FACSCalibur or a FACSCanto 
(BD), and data were acquired with CellQuest Pro or FACSDiva software, re- 
spectively, and analyzed with Flowjo software (Tree Star). 

Statistical analysis. Comparisons of phagocytic activity of macrophages of 
different genotypes were performed using one-way ANOVA and Bonferroni's 
multiple comparison post-test (a level 0.05) or two-tailed unpaired Student's 
t test using Prism software (GraphPad Software) . The statistical test, p-values, 
and the n for each statistical analysis are indicated in each of the correspond- 
ing figure legends. 

Online supplemental material. Figs. SI and S2 illustrate the genomic 
characterization of experimental mice. Tables SI and S2, included as separate 
PDF files, list relevant genes identified by RNA sequencing. Table S3 lists 
allelotypes of experimental mice. Table S4 lists primer sequences and restric- 
tion enzymes used for genetic analysis. Online supplemental material is avail- 
able at http://www.jem.org/cgi/content/full/jem.20131274/DCl. 
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